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We examine the spectrum of inflaton fluctuations resulting from any given long period of expo- 
nential inflation. Infrared and ultraviolet divergences in the inflaton dispersion summed over all 
modes do not appear in our approach. We show how the scale-invariance of the perturbation spec- 
trum arises. We also examine the spectrum of scalar perturbations of the metric that are created 
by the inflaton fluctuations that have left the Hubble sphere during inflation and the spectrum of 
density perturbations that they produce at reentry after inflation has ended. When the inflaton 
dispersion spectrum is renormalized during the expansion, we show (for the case of the quadratic 
inflaton potential) that the density perturbation spectrum approaches a mass-independent limit as 
the inflaton mass approaches zero, and remains near that limiting value for masses less than about 
1/4 of the inflationary Hubble constant. We show that this limiting behavior does not occur if one 
only makes the Minkowski space subtraction, without the further adiabatic subtractions that involve 
time derivatives of the expansion scale factor a{t). We also find a parametrized expression for the 
energy density produced by the change in a{t) as inflation ends. If the end of inflation were suffi- 
ciently abrupt, then the temperature corresponding to this energy density could be very significant. 
We also show that fiuctuations of the inflaton fleld that are present before inflation starts are not 
dissipated during inflation, and could have a significant observational effect today. The mechanism 
for this is caused by the initial fluctuations through stimulated emission from the vacuum. 
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I. INTRODUCTION 

Cosmological inflation predicts an amplification of 
quantum fluctuations that perturbs the background ho- 
mogeneity of the universe and leads to regions of space- 
time of nonuniform density |l|. The dispersion of the 
quantized inflaton fluctuation field, 54), is taken as a 
measure of the inhomogeneity of the field during infla- 
tion. This inhomogeneity perturbs the gravitational field 
during inflation, and its perturbations set the initial con- 
ditions for the acoustic oscillations of the plasma and 
matter that are present after the inflationary era. 

The dispersion of (50 has an ultraviolet divergence that 
results from integrating over all modes at a given time 
during inflation. As the highest frequency modes corre- 
spond to wavelengths that never exit the de Sitter hori- 
zon during inflation, one may think that such modes 
can be ignored or dealt with by means of a cut-off or by 
standard curved-spacetime regularization and renormal- 
ization methods without observable consequences. How- 
ever, 0] showed that when the method of adiabatic regu- 
larization is used to renormalize the dispersion of 5(j), the 
dispersion spectrum at wavelengths that have left the do 
Sitter horizon is significantly affected. Adiabatic regu- 
larization in Robertson- Walker universes has been shown 
to give the same result as other forms of regulariza- 
tion and renormalization, including point-splitting regu- 
larization; and the adiabatic condition has been shown to 
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be closely related to the Hadamard condition in curved 
spacetime. Because of the time-translation properties of 
the universe during inflation, renormalization of the spec- 
trum of inflaton fluctuations does not alter the near scale- 
invariance of the spectrum of perturbations as they en- 
ter the Hubble horizon of the post-inflationary universe. 
However, renormalization may alter the relation between 
the magnitude of the inflaton perturbation spectrum and 
the implied value of H during inflation. This would be of 
importance in testing theories that predict the value of 
H directly. The effects studied in [2] have been further 
elucidated in [1, @ . 

Is there a way to obtain the dispersion spectrum of 
inflaton fluctuations by using only the well-tested tech- 
niques of quantum field theory in Minkowski space with- 
out appealing to renormalization in curved spacetime? 
Such a method was developed by one of the authors (LP) 
in the early 1960's in the first part of his Ph. D. thesis 0]. 
First we explain how and why the method works. Then 
we trace the evolution of the quantized inflaton fluctua- 
tion field on the background metric, without considering 
the metric perturbations and the density perturbations 
that they lead to after reheating. In that case, we can use 
the method to obtain the dispersion spectrum and to ex- 
amine the scale invariance of the spectrum resulting from 
a long period of inflation. We also demonstrate how the 
method avoids both infrared and ultraviolet divergences 
in the dispersion of the inflaton flucutation fleld. 

Then we turn to processes involving the production of 
scalar perturbations of the metric by the modes of the 
inflaton fluctuation field after they exit the infiationary 
Hubble sphere; and to the spectrum of initial density 
perturbation they create upon reentry into the Hubble 
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sphere after the end of inflation. For these considerations 
we consider the consequences of using only the Minkowski 
space regularization, as compared with regularizing in- 
finities of the inflaton dispersion that involve derivatives 
of a(t) duriing the expansion, and using the regularized 
inflaton dispersion at the time when they induce signifi- 
cant metric perturbation after leaving the inflating Hub- 
ble horizon. We find that there are significant differences 
in the magnitudes of the density perturbations and their 
dependence on the effective mass of the inflaton fluctu- 
ation field. In particular, using the regularized disper- 
sion leads to initial density perturbations that approach 
a nonzero value in the limit of zero inflaton mass, and 
that are almost independent of the inflaton mass if that 
mass is less than a significant fraction of the inflationary 
Hubble constant. 

In 0, LP showed that an expanding universe creates 
elementary particles. He assumed that the quantized 
field is evolved by the generally covariant field equation 
in the expanding universe. If the quantized particle field 
was expanded in terms of mode function solutions of the 
field equation and the creation and annihilation opera- 
tors were defined as usual in terms of the coefficients of 
the mode functions, then he found that the particle num- 
ber density created in a given mode was finite, but when 
summed over all modes the particle number density had 
an ultraviolet divergence. This raised the questions, (1) 
should the particle number operator that yields the par- 
ticle number per unit physical volume be renormalized 
during the expansion of the universe, and (2) is it possi- 
ble to prove, using only the known and tested quantum 
field theory in Minkowski space, that the density of par- 
ticles created during the expansion of the universe, as 
predicted by quantum field theory, is actually finite? 

To answer the latter question, he considered a gen- 
eral expansion^ of the universe that started smoothly 
from Minkowski space in the distant past and after ex- 
panding in an arbitrary smooth manner approached a 
Minkowski space again in the distant future. It is rea- 
sonable to assume that the density of real particles cre- 
ated during the expansion is not greatly disturbed by the 
gradual and smooth joining to the early- and late-time 
Minkowski spaces. He showed that for such asymptoti- 
cally Minkowskian expansions of the universe, the num- 
ber density of particles created by the expansion of the 
universe was finite. No renormalization of the particle 
number in the curved spacetime during the expansion of 
the universe was used to obtain this result. Only the stan- 
dard definition of particle number in the initial and final 
Minkowski spaces was used. He employed general math- 
ematical theorems to prove that the particle number in 
a given co-expanding volume is an adiabatic invariant, 
and that the total number of particles created in the co- 
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expanding volume by the expansion of the universe from 
the initial to the final Minkowski space is finite when 
summed over all modes. This was proved by considering 
the Bogoliubov transformation, aj: = akAj:+f}lA^ that 
relates the annihilation operator for particles in mode 
k at early times to the annihilation and creation oper- 
ators at late times. He showed that the quantity |/3fcp 
that determines the average number of particles, created 
in mode k from the early-time vacuum by the expansion 
of the universe, vanishes faster than any inverse power of 
k, as k oo. Consequently, the average number of cre- 
ated particles, summed over all modes, is finite as mea- 
sured in the late-time Minkowski space. This result is 
independent of renormalization in curved spacetime. 

He addressed the former question (about renormaliza- 
tion of the particle number operator during the expan- 
sion) by considering the properties of a device that mea- 
sures the number of particles per unit physical volume. 
As shown in detail in [y| and summarized in 7] , the natu- 
ral assumption that the measuring instrument is not able 
to measure a particle number that has very fast small os- 
cillations, together with the requirement that the number 
measured in a given physical volume should be an inte- 
ger, leads to a renormalized definition of the Hermitian 
number operator corresponding to the quantity actually 
measured by such an instrument. This physically rele- 
vant renormalized number operator in an expanding uni- 
verse is the "bare" number operator obtained from the 
mode function expansion of the field, renormalized by 
making the adiabatic subtractions to second adiabatic or- 
der. The method of adiabatic regularization was further 
developed and applied to the energy-momentum tensor 
with Fulhng and Hu in f^. 

In order to obtain unambiguous results, independent of 
regularization and renormalization in curved spacetime, 
for the dispersion of the inflaton field resulting from a 
long period of inflation, we will consider a model of the in- 
flationary universe that starts smoothly from Minkowski 
space in the distant past and, after undergoing any given 
number of e-foldings of inflation, approaches a Minkowski 
space again in the distant future. (We do not take up per- 
turbations of the background metric that are produced 
by the created inflaton fluctuations until later in this pa- 
per.) We first take the initial state of the quantized infla- 
ton perturbation field to be the Minkowski vacuum. The 
asymptotic Minkowski space in the distant past may be 
regarded as a convenient way of specifying unambiguous 
initial conditions that can be interpreted without refer- 
ence to curved spacetime. The asymptotic Minkowski 
space at late times plays a similar role in permitting 
us to unambiguously interpret the spectrum of inflaton 
perturbations at late times after any given number of 
e-foldings of inflation. As explained above in our discus- 
sion of [y, 0], we must do the joining at early and late 
times as smoothly as possible to obtain results that do 
not have an ultraviolet divergence when summed over all 
modes. Our present method, in which a{t) is asymptoti- 
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cally flat at early times, also avoids infrared divergences, 
such as those that were found in the treatment of gravi- 
ton production from inflation by [9|. The infrared and 
ultraviolet divergences found by [lO| are also absent. 

In Section ini we specify the inflationary spacetime by 
joining together a composite scale factor built of seg- 
ments for which the solution to the evolution equation of 
the inflaton fluctuation field is known analytically. Be- 
cause observation [ll|, [l^] favors a period of nearly ex- 
ponential inflation with a slowly changing inflationary 
Hubble constant, iJinfl, we include as the middle part 
of our composite scale factor a region of exponential in- 
flation. We take Hi^a as a constant, but we can incor- 
porate a slow change by using the usual adiabatically 
adjusted solution in which the mode-function solutions 
are adjusted by neglecting the time-derivatives, but in- 
cluding the slow change of Hind- The initial and final 
asymptotically Minkowskian segments have a number of 
adjustable parameters that allow us to choose any value 
of ffinfl and any number of e-folds of inflation, while join- 
ing the scale factor a{t) continuously and with continuous 
first and second derivatives. That is, a{t) will be a 
function. This is the minimum degree of continuity of 
a{t) that, in general, will give a finite energy density in 
the late-time Minkowski space. We can also adjust pa- 
rameters that determine how quickly inflation ends. In 
Section IIIIl we obtain the asymptotic conditions on the 
modes of the inflaton perturbation field. 

In Sect ion [TVl we discuss the average number of quanta 
of the inflaton perturbation field for a pure state (the 
early-time vacuum) and for a statistical mixture of states 
having different numbers of particles present at early 
times. We discuss a surprising feature that results from 
the stimulated creation of particles if there are particles 
present at early times. This effect remains significant at 
late times, particularly at scales that may be relevant to 
the large scale structure of the universe. Without amplifi- 
cation by stimulated emission, the initial inflaton fluctua- 
tions (i.e., particles or quanta of the inflaton fluctuation 
field) would be dispersed by the inflationary expansion 
and would have negligible effect. We will take this up 
further in a later paper. In standard treatments of in- 
fiation, it is assumed that the initial conditions have no 
significant effect after a sufficient number of e- foldings. 
This statement, although it seems intuitively obvious, is 
not correct unless the initial state of the inflaton fleld is 
the vacuum state. 

In Section |Vl we solve the evolution equation for quan- 
tized fluctuations of the inflaton field for the above class 
of scale factors a{t). We do this by matching the analytic 
solutions for the modes of the inflaton field perturbations, 
(5(/>, and their time-derivatives at the joining points where 
we spliced together the different segments of the scale fac- 
tor. We introduce the general evolution equation for Scj) 
with a constant effective mass and focus initially on the 
exact solutions of the massless case. In Section I VII we 
find the number of particles created in each mode for in- 
fiaton fluctuations of effective mass, and we discuss the 



effect of discontinuities in a(t) and its derivatives on the 
average number of inflaton perturbation quanta created 
in each mode. 

In Section IVIIl we obtain the dispersion spectrum of 
the massless inflaton fluctuation field. As mentioned 
above, there are no infrared or ultraviolet divergences 
in these quantities. We also find effects that are evident 
in our dispersion spectrum and depend on the phases of 
the inflaton perturbation modes at the time that inflation 
begins in our model. 

In Sections IVIIII and IIXI we find the corresponding 
results for quantized inflaton fluctuations of non-zero ef- 
fective mass. In Section |Xl we discuss the dependence 
of the spectral index on the effective mass of 6(j). We 
also investigate the effect of the total duration of infla- 
tion on the scale-invatiance of very long wavelength per- 
turbation modes. In Section IXIl we discuss the initial 
density perturbations we would expect to be produced 
after reheating, and we show that there is a significant 
difference when we employ renormalization of the infla- 
ton dispersion during the expansion, as compared with 
employing only the Minkowski space renormalization. In 
Section IXIIi we use our parametrized scale factor to cal- 
culate the contribution of the changing gravitational field 
to reheating (see also [Tsj). 



II. COMPOSITE SCALE FACTOR 

We consider the metric 

ds^ = dt^ - a^{t){{dxf -f {dyf + {dzf). (1) 

The time t will run continuously from — oo to oo. The 
scale factor a{t) will be composed of three segments. 
Our scale factor will generally be C^, i.e., a continuous 
function with continuous first and second derivatives ev- 
erywhere, including at the joining points between seg- 
ments.^ The initial and final segments are asymptoti- 
cally Minkowskian in the distant past and future, respec- 
tively. The middle segment is an exponential expansion. 
We choose specific forms for a{t) in these segments that 
have exact solutions of the evolution equations for infla- 
ton quantum fluctuations of zero effective mass. 

We emphasize that the initial and final asymptotically 
fiat regions permit us to unambiguously interpret our re- 
sults for free fields without having to perform any renor- 
malization in curved spacetime. The final asymptotically 
flat region will not signiflcantly affect the result obtained 
for the spectrum of inflaton perturbations created by the 
inflationary segment of the expansion. The initial asymp- 
totically flat region should have a negligible effect on the 



^ We also briefly consider scale factors that are only or C" at 
the joining points. 
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spectrum resulting from a long period of inflation.^ If 
there are cases in which no inflaton perturbations are cre- 
ated by the period of exponential inflation, then the in- 
flaton perturbations would result from the initial asymp- 
totically flat region, possibly amplified by the long period 
of inflation . 
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FIG. 1: Scale factor, a(t(r)), and dimensionless Hubble pa- 
rameter, sH{t{T)) — sa~^da/dt = sa~'^da/dT, of Eq. ((2)1 with 
ai = 1, a2 — 2, b = 0, and s = 1. Note in the graph that 
the maximum of H occurs at a value of a(t(r)) closer to ai 
than to a2. In both the case where 02 3> ai and the case 
where 02 — ai, -ffmax occurs at a value of the scale factor 
where a(f(r)) ~ ai. The relationship between proper time 
and r-time is given by Eq. (|3}. 

We base each asymptotic segment on a scale factor of 
the form, 

a(t(r)) = 

1 

|a/ + e^/^[(a24 _ a,^){e^/^ + 1) + b]{e^/^ + Ij-^j ' ,(2) 

where r is related to the proper time as defined through 
Eq. ©, by 

dr = a{ty^dt. (3) 

The form of the scale factor in Eq. ^ is based on 
the form of the index of refraction used by Epstein to 
model the scattering of radio waves in the upper atmo- 
sphere and by Eckart to model the potential energy in 
one-dimensional scattering in quantum mechanics [14| . 
It was first used in the cosmological context by Parker 
(m to model a{t). As can be seen from Fig. [TJ this scale 
factor approaches the constant ai at early times and the 
constant a2 at late times, and the constant s determines 
roughly the interval of r-time for a{t) to go from ai to 



^ We do find remnants of the early initial conditions in the late- 
time inflaton dispersion spectrum, which we discuss later. 
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FIG. 2: Scale factor plotted versus dimensionless time. This 
illustrative example summarizes our notation using a moder- 
ate expansion of ~ 2 e-folds. The scale factor, a{t), is con- 
tinuous, as are d{t) and a{t). In this case, the parameters for 
the initial asymptotically flat segment are an — 1, a2i = 2, 
and Si — 1. The free parameters of the final asymptotically 
flat segment are 02/ = 9 and ai/ —6. The asymptotically 
flat scale factor of the initial region joins the exponentially 
expanding scale factor of the middle region at a time ti in 
t-time and Ti in r-time. The exponentially expanding scale 
factor of the middle region joins the asymptotically flat scale 
factor of the final region at a time t2 in t-time and r| in r'- 
time of the final segment, where a prime is used to distinguish 
between the r-times of the initial and final segments. 

02. (A sufficiently large magnitude of b would produce a 
bump or valley in a{t).) The parameters oi, 02, 6, and 
s are different in the initial and final asymptotically flat 
segments. Where confusion would arise we will include 
subscripts i in the initial set of parameters and / in the 
final set of parameters. The equation for the middle (in- 
flationary) segment of our composite scale factor is given 
in terms of proper time by 

a(i) = a(ti)e^'""(*-*i), (4) 

where Hiad is the constant value of H{t) = a^^da/dt 
during the exponential expansion of the middle segment. 

We define the quantity = In {a2f/aii)- When there 
is a long period of exponential growth, is essentially 
the number of e-foldings of inflation. Typically, Ne will 
be about 60. Within the final asymptotically flat scale 
factor, the ratio of 02/ to aif determines how gradually 
the exponential expansion transitions to the asymptoti- 
cally flat late-time region. (For example, this ratio might 
be I e-fold, which we would consider to be relatively grad- 
ual, or it might be 1.0001, which we would consider to 
be relatively abrupt.) 

With our choices of a{t) in the three segments, we are 
able to join them so that a{t) and its first and second 
derivatives with respect to time are continuous every- 
where. This requires that we join the exponentially ex- 
panding segment, in which H{t) has the constant value 
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ffinfl, to the initial and final segments at the times when 
H{t) is a maximum. This maximum value must equal 
ifinfl- A simple power law form of the scale factor, such as 
that of a radiation-dominated universe, could not be used 
to simultaneously maintain the continuity of the scale 
factor and its first and second derivatives when matched 
directly to the inflationary segment of exponential expan- 
sion. 

With bi = and 6/ = 0, we then find the following 
expressions. The time at which the first segment joins 
to the exponential segment is 

C, 



Si In ' " 



The constant a{ti) in Eq. ^ is 



aih) 



1/4 



(5) 



(6) 



Because the maximum value of H(t) in the first segment 
must equal i?infl, we find that 



\l/4 



where 



C,= 



(7) 



(8) 



Once we choose values for ai/ and 02/, the remaining 
constants are determined to have the following values: 



Sf = 



23/4 (_«^; 



^2/ 



llaj^^ -3a,f* + Cf) H, 



Mnfl 



X (-3ai/ - 302/ -I- C/) 



1/4 



X (3ai/ - 3^2/ + Cf) . (9) 

Denote the parameter r of Eq. ([2]) as r' in the final seg- 
ment. At the time t'j when the exponential segment joins 
to the final segment, we find that 



t'j = Sf In 



(10) 



^2/ 



The corresponding proper time t at which the exponen- 
tial segment joins to the final segment is 



t2 = 



1 



■In 



— 3a]^j4 



302/ + Cf 



where 



infl 



Cf = 



3a2»'' 



+ ti 



(11) 



(12) 



See Fig. [5] for a schematic diagram of how we match 
our segments of the scale factor together. 



III. EARLY AND LATE ASYMPTOTIC 
CONDITIONS ON 5<p 

Consider an inflaton field composed of a spatially ho- 
mogeneous term plus a first order perturbation, 



(13) 



We investigate, in units of ft = c = 1, a minimally- 
coupled scalar field that obeys the evolution equation: 

3 

dt^S(P + 3HdtS(p - a-'^{t) d,^6(f) + m(?!)(°))250 = 0. 

1=1 

(14) 

The mass term is related to the infiationary potential by 



rf((/.(0))2 



(15) 



For simplicity, we take 771(0^°-*)^ as a constant, m^. 

The quantized field S4> can be written in terms of the 
early time creation and annihilation operators, Al and 
A^, as 

k k 

where^ 

/- = y-le^fc-->fc(t(T)). 
The function ^->k{t) satisfies 



(16) 



(17) 



d''Mt) + ^HdtMi) + ^7-rV'fe(t) +"iVfc(t) = 0, (18) 
a (*) 

where k = 27rn/L, with n an integer. Because the cre- 
ation and annihilation operators in Eq. (1161) correspond 
to particles at early times, we require that i/jfe satisfies 
the early-time positive frequency condition 



lim ?/;fe(t(T)) 



T^ — OO 



yj2aliUii{k) 



(19) 



where LOii{k) = \J {k/au)"^ + m?. 

At late times, this solution will have the asymptotic 
form 



lim Mi{T')) 



^2al^uj2f{k) 
-|-/3fce*°2/"2/(fe)T-' 



(20) 



4 For simplicity, we are imposing periodic boundary conditions 
upon a cubic coordinate volume, V = . In the continuum 
limit L would go to infinity. 
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where uj2f{k) = y/ {k/a2fY + 'm-'^- Here, one has 
|«fc|'-|/3fc|' = l, 



(21) 



from the conserved Wronskian of Eq. (jl8[) . The quantity 
|/3fc|^ is the average number of particles in mode-fc cre- 
ated by the expansion of the scale factor from a state that 
initially has no particles [1, 0] ■ Because our scale factor 
is asymptotically Minkowskian, the meaning of particles 
at early and late times has no ambiguity. The late-time 
creation and annihilation operators, at and ag, are re- 
lated to the early-time creation and annihilation opera- 
tors through a Bogoliubov transformation: 



(22) 



IV. 



PARTICLE NUMBER IN PURE AND 
MIXED STATES 



As noted above, if no particles are present at early 
times, then the average particle number at late times in 
mode k is 



= (0|aU,,|0) = |/3fcf 



(23) 



where |0) is the state annihilated by the early-time anni- 
hilation operators A^. 



Let us define ^^^"'^^ 


2 


(%)2 0)u„ - \fj:\\ Here, 


"un" refers to unrenormalized values. In the continuum 


limit, this reduces to 




= (27r)-3|V.fc|2. 



Later, we will find that renormalization is necessary 
and may have a significant effect on the magnitude of 
the dispersion, even for modes that leave imprints on 
the CMB and large-scale structure that can be observed 



in the present universe, 
where 



Let 



(0 {5h? 0) 



^ (0 |(%)2 Ore ^ (0 |(J0(X))2| Ore. (24) 



The value of (0 |((50(x))^| 0)ic would diverge without 
renormalization. In Minkowski space, the renormaliza- 
tion would be equivalent to subtracting the vacuum zero- 
point contributions. In the expanding universe there 
are also divergent contributions coming from the time- 
dependence of the scale factor a(t), and to take those 
into account we must use a curved spacetime renormal- 
ization method such as adiabatic regularization or point- 
splitting Hadamard regularization. We will take this up 
again later in Section Ixj 

If the initial state of the universe were not the vacuum, 
but instead were a statistical mixture of pure states, each 
of which contains a definite number of particles at early 
times, then we would find the analog of stimulated emis- 
sion, where the initial presence of scalar particles tends 



to increase the number of scalar particles created by the 
expansion of the universe 0, 0] : 



(25) 



Here, {Nj:)^^^ is the average particle number in mode k 
at late times, and {N^^) is the average number of parti- 
cles in mode k at early times. Because of this stimulated- 
emission effect, the initial presence of particles could lead 
to larger inflaton perturbations than would be the case 
with an initial Minkowski vacuum state. Furthermore, 
it could affect the scale invariance of the inflaton pertur- 
bation spectrum at late times. The particles that were 
present at early times disperse so that their relic density 
becomes negligible as a result of inflation. This is evident 
from the first term on the right of Eq. (|25|1 . However, the 
density of quantized inflaton perturbations created as a 
result of the inflationary expansion remains significant 
and is determined by the second term. In that term, the 
factor of |/3fc|^ multiplying {Nj^°) makes the effect of the 
initial number of particles present in the coordinate vol- 
ume significant if (^g") is larger than or of order 1. 
This effect could conceivably lead to observable conse- 
quences for the large scale structure of the universe, and 
will be considered in later work. In the rest of this paper, 
we assume our initial state is asymptotically a Minkowski 
vacuum. 



V. JOINING CONDITIONS FOR Vfc 

Consider a spacetime composed of three segments of 
the scale factor, a{t), in a homogeneous background met- 
ric given by Eq. ([T]). For an example, see Figs. [2] and [3l 
The first and second segments are joined at the time ti , 
and the second and third segments are joined at the time 

t2- 

We have two linearly independent solutions to the evo- 
lution equation in both the second segment, with solu- 
tions hi{t) and h2(t); and the third segment, with solu- 
tions gi (t) and g2 (t) ; for a total of four separate functions. 
These functions are multiplied by constant coefficients 
that we must determine. During the second segment, 
from ti to t2, we have: 

Mt) = Ahi{t)+Bh2it), (26) 
^',{t)^Ah[{t)+Bh'2{t). 

For t > t2, we have: 

Mt)=C9iit)+Dg2{t), (27) 
^Piit) = Cg[{t) + Dg'2{t). 

If we require that ipkit) and tp'kit) be continuous at ti 
and t2. This imposes 4 matching conditions: 

Ahiiti) + Bh2ih)^Mh), (28) 
Ah[it,) + Bh'^{t,) ^ ij^ih), 
Cgi{t2) + Dg2{t2) - Ahi{t2) + 5/13(^2), 
Cg[{t2) + Dg'^{t2) = Ah\{t2) + Bh'^{t2). 
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The calculation given in Appendix 1X1 then shows us that 
1 



C = 



{a'i92 - 9192) 

fcl"2 



+ 



hlh'2 

iJjkiK 



{h'i92 



h^hi 



h2h[ 



{h'292 



h292)t=t2 



t=ti 



(29) 



D 



(5251 - 929'i)t=t2 

h[h2 — hih'2 
■0^1 /ii - ipkih'i 



h^hi — h2h[ 



{h'i9i 
ih'291 



hi9'i)t=t2 
h29i)t=u 



(30) 



where ipki = V'fc(^i) a-nd V'fei = '0fc(^i)- We find if^ki and 
'0^]^ from the solution to the evolution equation in the 
initial asymptotically flat segment of the scale factor. In 
the massless case, this solution is given by Eq. ([5^ below. 
The functions hi{t) and h2{t) are to be related to the 
evolution equation solutions in the inflationary middle 
segment of the scale factor, which are given in Eqs. pS]) 
and ([53]) below. In terms of those solutions, one finds 
that A and B in Eq. are given hy A = E{k) and 
B = F{k). Similarly, the functions gi{t) and 32 (i) are 
related to the solution of the evolution equation in the 
final asymptotically flat segment of the scale factor. The 
latter solution is given in Eq. (|34p . which allows us to 
specify that C = A^i(fc) and D = N2{k). 



VI. MASSLESS PARTICLE PRODUCTION 

We will first consider the case, ra — Q. Rewriting the 
evolution equation, Eq. p8)) . in terms of r instead of t 
leads to 



(31) 



For the first segment of our composite scale factor, the 
solution of ([?T|) having positive frequency form jigp at 
early times is the hypergeometric function (isl . ITT 



Mtir)) 



1 



, =e ^'"'^'^'^ F(-ikau^Si + ika2i'^Si, 

—ikaii^Si — ika2i^Si] 1 — 2ikaii'^Si] — e ), 

(32) 



where F(a,b;c;d) is the hypergeometric function as de- 
fined in [13, see 15.1.1]. 



For the exponentially expanding segment of the scale 
factor in the massless case 



F{k)H 



(2) 



E(k)H^^^ 

2 

k 



i{t)Hi 



infl 



(33) 



where i?^^^ and iJ^^^ are the Hankel functions of the first 
and second kind. The variables t and r are related by 
Eq. ([3]). The coeflicients E{k) and F{k) are determined 
by the matching conditions of the first joining point at 
t = ti. We note that the finite period of exponential 
inflation lacks the full symmetries of a de Sitter universe. 
In the pure de Sitter CBtSG, clS shown in [l^, the fc = 
mode has to be chosen in a special way to avoid infrared 
divergences. For our a{t), infrared divergences do not 
arise (see Sec. IVI|) . 

For the final segment of our composite scale factor, 
the solution of the evolution equation (1311) is a linear 
combination of hypergeometric functions jl5l.[l6j: 

Mt{r')) = Ni{k)e-''''''f'^' Fi-ikaif^sf + ika2fhf, 

—ikaif'^Sf — ika2f'^Sf \ 1 — 2ikaif^Sf] —e"f ) 
+N2ik)e''"''f^''' F{ikai/sf + ika2/sf, 

ikaif^s — ika2f'^Sf; 1 + 2ikaif^Sf; —e'l ), 

(34) 

where the coefficients Ni{k) and N2{k) are determined 
by the matching conditions of the second joining point 
at t = ^2- 

An example of the evolution for a particular mode 
is plotted for a specific choice of parameters using our 
composite scale factor in Fig. [31 We emphasize that our 
considerations below apply to any length of the middle 
exponentially inflating segment of the expansion of the 
universe, not just to the one used in Figs. [2] and [3l 

For our choice of the flnal asymptotically flat segment 
given by Eq. where we use Eq. ([M)) to define our 
functions gi(t) and 52 (i) in terms of the relationship 
V-'fe(i) — Ni9i{t{T)) + N2g2{t{T)), we find the coeflicients 
Q!fe and f3k of E g. (120| from the large argument asymp- 
totic forms [Tsl. [igI. Il7| . With bf = 0, ci = iksfa^^, and 
C2 = iksfa2f, we have 



ak 



D r(i 



c r(i - 2ci)r(-2c2) 
r(i - ci - C2)r(-ci - C2) 

2ci)r(-2c2) 



and 



Pk 



r(i + ci- C2)r(ci - C2) 



cr(i-2ci)r(2c2) 



(35) 



^kcii2 j: 



r(i-ci + c2)r(-ci + c2) 
D r(i + 2ci)r(2c2) 



r(i + ci + c2)r(ci + c2) 



(36) 
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FIG. 3: (color online). A solution to the evolution equation 
for a particular Fourier mode ipk{t) is plotted versus dimen- 
sionless time for the same composite scale factor a{t) shown 
in Fig. [21 The real and imaginary parts of the dimension- 
less quantity y/^xjjkit), and its magnitude are all plotted for 
fc = 2. 



Recall that C and D and the functions gi{t) and g2{t) 
were defined in Sec. |Vl As the Wronskian of the evo- 
lution equation is conserved with our joining conditions, 
Eq. (PT|) should hold to arbitrary accuracy. As a useful 
check on our calculations, we verified that in all the cases 
considered in this paper, the relation |afc|^ — |/3fc|^ — 1 is 
satisfied to at least 500 significant figures. 

The quantity |/?fcp is the average number of particles 
present in mode k at late times. These particles are cre- 
ated by the expansion of the universe through N^. e-folds 
from a state having no particles present at early times. 

We use the dimensionless variable 



92 



(37) 



where k is the wave number, 02/ is the constant scale 
factor that is approached at very late times, and H-i^a is 
the constant value of {d{t)/a{t)) during the exponential 
expansion of the middle segment. It is convenient to 
express our results in terms of the dimensionless quantity 
92- For example, for Ne — 60, as in Fig. IH the graph is 
the same if we change the range of k, and the values of 
02/ and iJinfl, while keeping the range of 52 unchanged. 
We will often write IPg^ |^ for the quantity |/3fc|^ with k = 

We define three regions of q2. Values of q2 < exp(— A^e) 
are in the "small-(72 region." Values of exp(—A''e) < 92 ^ 1 
are in the "intermediate-g2 region." Values of 1 < 92 are 
in the "large-g2 region." 

When a(t) is at least C^, i.e. when HinU is contin- 
uous, we find numerically that the particle production 
per mode in the small-g2 region, (52 ^ e~^=), is to good 
approximation given by 
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FIG. 4: Average late time particle number per mode, (|/3ij2 P), 
versus 52 = k/{a2fHina) for 60 e-folds of inflation. Two cases 
are plotted for the massless case based on the behavior at 
the matching conditions: the scale factor continuous in 0th, 
1st, and 2nd derivatives (C^); and the scale factor contin- 
uous in 0th and 1st derivatives (C^). Note that in the 
case, 1/3,2 1^ transitions from a ^2"^ dependence at the end of 
the intermediate-g2 region all the way to a ^2"^ dependence, 
temporarily parallel to the large-g2 regime, before settling 
down into its ultraviolet ^2"* behavior. For the wiggles near 
the transition from the small-g2 region to the intermediate-g2 
region at ^2 = e~^^ , compare with the graph of the dispersion 
spectrum in Fig. [S] 



We also find this to be the case, analytically, by taking 
the limit fc — s- 0. For at least a moderate number of 
e-folds, this simplifies to 



1/3, 



92 1 



1 



.2N, 



(39) 



The dependence in the intermediate-(72 region (e ^= 
92 ^ 1) for the or massless case is 



< 



1 



(40) 



/3,2 = sinh[A^e 



(38) 



When Nf. is finite, with our composite scale factor there 
are no infrared divergences. For infinite inflation, where 
Nf, — > 00, we find the infrared divergences of a de Sitter 
universe. This problem is resolved for a true de Sitter 
universe in [isj . Our composite scale factor is different 
from a purely de Sitter universe in that our initial condi- 
tions are specified by our initial asymptotically flat region 
of the scale factor. 

Discontinuities in the derivatives of the scale factor at 
the matching points cause increased particle production 
for large-g2, (i.e., 92 This is evident in Fig. [H As 
shown by Parker in [g, |7|] , the particle number created by 
the expansion of the universe is related to an adiabatic 
invariant of the harmonic oscillator, and for such an oscil- 
lator the continuity of the frequency and its derivatives 
is related to the change in the adiabatic invariant p^ . 
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The increased particle production associated with dis- 
continuities in the derivatives of the scale factor a{t) is 
a consequence of this relation. This result is also seen in 

For the case, where the scale factor and H = 
d(t)/a{t) are both continuous, the large-(72 region goes 
like 

= (41) 
For the case, where the scale factor and H — d{t)/a{t) 
and H{t) are all continuous, the large-(72 region goes like 



"692 



(42) 



Here and ng are constant coefficients, with n4 ~ ng ~ 
0(1) for a gradual end to inflation. For a sufficiently 
abrupt end to inflation, 77.4 and ng can be made to be 
arbitrarily large. See Sec. I XIII 

When H(t) is not continuous, we find quite a differ- 
ent behavior in the C° case. The evolution equation, 
Eq. (HH), may be written 0] 



ait) 



3d{t) 
2 a(i) 



Mt) = 0. 

(43) 

At the discontinuity in d(t) if we express the jump as 
a step function, then the form of d{t) picks up a delta- 
function contribution. Thus, there is a finite jump in 
dipkit)/dt across the discontinuity. In the C° case, 1/3,2 P 
is proportional to 93"^ in the small- and large-g2 regions, 
and it is proportional to 92"^ '^^ the intermediate-g'2 re- 
gion. A C° scenario would suffer from these problems in 
addition to the divergences mentioned earlier, hence we 
will not consider it further. 

For a non-composite scale factor composed of one 
asymptotically flat scale factor defined by Eq. for 
large values of q2 the value of |/3g2p falls off faster than 
any power of q2, and in terms of k we have: [l^ [l^ 

sin^ - VTTIk^]) + sinh2[7rfcs(ai2 - a^^)] 

(44) 

In the limit that A: — > for the case of the scale fac- 
tor of Eq. which is asymptotically flat at early 
and late times and has no exponential segment, we flnd 
that limfe^o |/3fe|^ = sinh^[iVe], where in this case is 
In (02/01). This is the same small-g2 limit for the average 
number of particles created per mode as we found above 
in Eq. 



VII. DISPERSION SPECTRUM 

The unrenormalized dispersion spectrum is (l6l . [21I 



^(fc/a2/)2 +m2 



, (45) 



where the expectation value is with respect to the state 
|0) having no particles at early times. 



0.015 




FIG. 5: (color online). Dispersion ( | S(l}'^ \ ) /^^infl given by 
Eq. (|50|l for our composite scale factor continuous in a(t), 
a{t), and a{t) over an expansion of 60 e-folds. The y-axis, 
( 1 ^'P^ I )q2 /-^infl) is shown multiplied by a factor of e'^"; 
and the x-axis, q2, is shown multiplied by a factor of e^" . 
When using this scaling, for a given set of values of an and a2i 
the region plotted in this graph would look identical for any 
number of e-folds larger than about 10. In the case of a2i = 
an + dan, where dan = 10~^^aii, we see marked peaks in the 
dispersion spectrum. When we change the parameters in the 
initial asymptotically flat region to a2i — lOaii, these peaks 
are somewhat damped, as shown. The ending conditions of 
the final asymptotically flat segment do not affect these peaks. 

We will flrst consider the massless case, m = 0. See 
Sec. [lX]for the massive case. If |/3fc|^ were in Eq. ([iSjl . 
then one would be left with the contribution of the late- 
time Minkowski space vacuum. As is standard in flat 
spacetime quantum field theory, we subtract off this con- 
tribution of the Minkowski vacuum, leaving the physi- 
cally relevant, renormalized dispersion. 



1 

1 



E 



^2/ 



i3 



E 



(46) 



where the subscript "re" refers to renormalized. In the 
continuum limit, this becomes 



1 



a2/(27r)3 
or with spherical symmetry, 

-2^ 

With k = q2a2fHinfi and dk = dq2a2fHiafi, we have 



(47) 



27r2 



q2\Pq2\ dq2. 



(48) 



(49) 
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FIG. 6: (color online). Dispersion spectrum Pi^Z-ffinfl given 
by Eq. (|52|l for our composite scale factor continuous in a{t), 
a{t), and a{t) over an expansion of 60 e-folds. The two cases 
considered are the same as explained in Fig. \5\ Similar peaks 
were found in a function related to the energy-momentum 
tensor in [2^ . see their Fig. 1]. 



reduces to this standard Minkowski space subtraction in 
flat spacetime. In Sec. IXI[ we consider the scalar per- 
turbations of the metric that are created by the infla- 
ton fluctuations. Those metric perturbations for a given 
mode are formed during inflation shortly after the infla- 
ton mode has exited the Hubble sphere. As explained 
in 0, [E] J we regard the metric perturbations as classical 
at the time when they are produced and assume that 
they respond to the renormalized dispersion of the in- 
flaton fluctuation field. This implies that the relevant 
adiabatic subtractions that influence the metric pertur- 
bation for a given mode are those evaluated near the time 
that the quantized inflaton fluctuations exit the Hubble 
sphere. 

The definition of the spectrum of inflaton perturba- 
tions, Vstj), given in [2^ is 



Jo 92 



(51) 



Therefore, we obtain the spectrum of inflaton fluctua- 
tions as 




27r2 



(52) 



20 40 
T r 60t 

Ln[q2e ] 

FIG. 7: (color online). Dispersion spectrum Vs^i/H^^f^ given 
by Eq. (|52|l for our composite scale factor continuous in a{t), 
d{t), and a{t) over an expansion of 60 e-folds. The cases 
considered are the same as those in Figs. [S] and [S] Here we 
plot the data from Fig. [6] on a In- In scale over a wider range 
of 52 for a gradual end to inflation, where ln(a2//ai/) ~ 1. 



The g2-component of the dispersion in the massless case 
is thus, 



We plot ( |(5</>2| )^JHI^^ in Fig. El and we plot 

Vs^/Hf^f^ in Fig. [H The peaks in these figures do not 
depend on the number of e-folds and are the result of the 
state |0) having no particles in the asymptotically flat 
part of the expansion at early times. In Fig. [71 we plot 
ln(P50/i7?jfj) versus ln(q2 exp(60)) for an expansion of 60 
e-folds. In this figure, we have plotted the full range of 
the inflationary part of the expansion. We see that in 
both of the cases plotted, where a{t), d{t), and d{t) are 
all continuous; and the case where only a{t) and a{t) are 
continuous; (0| Scf)'^ [0)^.^ is flnite without the need for any 
renormalization beyond subtracting off the Minkowski 
vacuum contribution at late times. If the first deriva- 
tive of the scale factor were not continuous, then the 
integrated dispersion would diverge. 



VIII. MASSIVE PARTICLE PRODUCTION 

In the case of a massive scalar field, the evolution equa- 
tion p8|) with a{t) given by the infiationary exponential 
of Eq. ^ , has the exact solution 



'92 



27r2 



(50) 



In the late-time Minkowski space, this subtraction is 
all that is necessary, but during the time when the uni- 
verse is expanding this subtraction alone would give an 
infinite dispersion when summed over all modes. Dur- 
ing the expansion, there are additional subtractions nec- 
essary. Adiabatic regularization [s^ includes those and 



V 4- 



+F(fc)ff^ 



where 



(53) 



(54) 
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The evolution equation for arbitrary a{t) can be writ- 
ten in terms of r (defined in Eq. ([3])) as 



2„6n 



(55) 



In the initial and final asymptotically flat segments of 
a(t) that are joined to the inflationary segment at early 
and late times, we do not have an exact solution of this 
equation for nonzero mass. Therefore, we used approxi- 
mations to obtain the plots in Fig. [8] and Fig. [9l In these 
graphs, we see that at low momentum there are many 
more particles with = 0.1 than with mu = 1 present 
at late times. 
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FIG. 8: (color online). Log-Log plot of the average number of 
created particles versus q2. The effect of three different masses 
is shown for an expansion of 60 e- folds. The beginning and 
end segments of a{t) are rather abrupt, with a2i = a\i(l + 
10"^^), a2f = aiie®", and ai/ = 0.9999a2/. Comparing the 
graph plotted here for mn ~ with the corresponding graph 
plotted for a gradual end to inflation in Fig. 3] one sees that 
the two graphs are the same for 52 ^ 1. In the present graph, 
however, there is a "stretched" region of 1/3,2 1^ oc g^'^ to the 
right of q2 — 1 that lasts until 52 — 10* before the ultraviolet 
behavior of l/Sqal^ °^ i^ seen. (For more discussion, see 
Sec. lXIll l This stretching is due to particle creation caused by 
the rapid change of a{t) from its value at the end of inflation 
to its final value of 02/- Both of the massive cases shown 
here produce more particles of low momentum than does the 
massless case. See also Fig. [5] 



IX. MASSIVE DISPERSION SPECTRUM 




-30 -20 -10 
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FIG. 9: (color online). The dependence of particle produc- 
tion (|/3q2p) on mass is shown for an expansion of 60 e-folds. 
This graph is different from Fig.|5]in that the transition from 
exponential expansion to the final asymptotic segment of the 
scale factor is more gradual, happening over about an e-fold. 
For values of 52 ^ 1, this graph is identical to that of Fig. [S] 
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FIG. 10: 



(color online). 
r2 



Comparison of dispersion spectrum, 
I )q2 /^infl) given by Eq. (|56|) and normalized to 1, for 
our composite scale factor continuous in a{t), a{t), and a(t) 
over an expansion of 60 e-folds for various masses. The values 
of ( I 34)^ I /^^infi were divided by the maximum value of 
the primary peak located at 52 — exp(— A^a) for each. To 
normalize these peaks, ( | Sfj)^ \ ) /H^^f^ was divided by the 



following factors 
2 



1.3 X 10 

for = 0.1, and 2.2 x 10" for = 1 



for the massless case, 2.3 x 10 

,4 r 2 



A calculation of the dispersion spectrum in the massive 
case leads to an equation analogous to Eq. (|50p : 



I2 



where 



12 



(56) 



(57) 



The dispersion spectrum is plotted for three different 
cases of rriH in Fig. [TOl We have found that, even in the 
massive case, the observed humps are dependent only 
upon the initial conditions. The shape of the curves 
is fixed above a moderate number of e-folds. We de- 
fine the variable J, such that the maximum value of 
( I 6(j)'^ I /H^nR fo'^ tl^s major peak, which is the peak 
located nearest to q2 = e~^', is J e^^"^)-'^^ in the mass- 
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FIG. 11: (color online). Dispersion spectrum Vs,i>/H^^fi given 
by Eq. (|52|l for our composite scale factor continuous in a{t), 
d{t), and a{t) over an expansion of 60 e- folds. The masses 
considered are the same as in Figs.|5]and[51 where the average 
number of created particles was plotted. 

less case and is J e'^'^^^^^^ in the massive case. Then, the 
normalization factor scales like e'^''^^^-'^^ in the massless 
case, as can be seen from Eq. (1501) : and the normalization 
factor scales like e^^~'^^-^'' in the massive case, as can be 
seen from Eq. (j56p . where we define the exponent P in 
the following way: 

\P,f^\q,-'' (58) 
in the region of intermediate-g2, (e'^' ^ 92 ^ 1), and 

IP.f^y''^ (59) 

in the small-q2 region, ((72 ^ e^^'). The exponent P is 
well described by a g2-independent value in the case of 
m = and in the case of 0.01 ^ mj^ <9/'i. 

In the massless case, P = 2, as can be seen by Eqs. ([39]) 
and (|40|) . The height of the major peak in the graph of 
the massless case in Fig. [TU] grows with an increasing 
number of e- folds as e^" , while the widths of the peaks 
narrow with an increasing number of e-folds as e"^". 
The area under an individual peak in the massless graph 
therefore does not change appreciably when changing the 
number of e-folds of expansion, provided there are at least 
a few e-folds of inflation. For the massive cases, we see 
that P = 2.93358 when to J = 0.1, and that P = 2.23607 
when tuh — 1. See Fig. [T^l 

We wish now to approximate the dependence of the 
configuration space dispersion (0| \0)j.^/H^^f^ on the 
number of e-folds. This dispersion is obtained from 
Eq. (|45p and its continuum limit. It is proportional to 
the area under the curves of ( | (50^ | /H^nR [IHl 
for the particular values of the mass shown. The main 
contribution to the area under each curve comes from the 
intermediate values of q2 (i.e., exp(— TVg) < 92 < !)■ For 



I ^ ^ ^ 1 

0.5 1 1.5 2 

FIG. 12: The dependence of the variable P, as defined in 
Eq. (|58|l . upon niH = m/Hinfi- The calculated data points 
shown lie on the curve P — \/9 — 4m ^ . Outside of the region 
plotted, however, P does not have a constant, g2-independent 
value. For rriH > 1.5, the argument, ^ (9/4) — m^, of the 
Hankel functions becomes imaginary, and l/S^jl^ oscillates 
with changing §2. For an example of a non-zero mass much 
smaller than Hi^a, see Fig. 1131 
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FIG. 13: Particle production as a function of q2 is plotted 
for 60 e-folds for both the massless case and the case of m = 
10~^''iifinfl, labeled as "m << H." As the mass decreases, the 
region of overlap of the solid curve with the dotted (m — 0) 
curve becomes larger and approaches the dotted curve in the 
limit as m ^ 0. The solid curve breaks away from the dotted 
curve when 52 < mH- The solid curve takes the constant 
value, (1/4)52^^", when 52 < mHexp(— A^e). In the region 
of ran exp{—Ne) < q2 < rnn, we see {k/a(t))^ ^ m? in the 
initial asymptotically fiat region and {k/a{t)Y ^ rn^ in the 
final asymptotically flat region. Between 52 — mHexp(— A'^e) 
and 52 — exp(— A'e), we see l/Jqsl'^ ^ 12~^^ ^^'^ between g2 — 
exp(— Ae) and §2 — rnn, we see l/Jqal^ ^ ^2~^ ■ ^ri light of 
these characteristics, a comparison of Eqs. (I50p and (I56|l can 
be made with consideration to where {k/a{t)Y 3> rn^ and to 
where {k/a(t)Y ^ nn? ■ Such an analysis shows that in the 
tiny mass limit of mn ^ 1, the dispersion spectrum reduces 
to the massless dispersion spectrum. 
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TABLE I: Configuration Space Dispersion 







0.01 <mH< 1.49 




( 1 1 l-P JP-3)N^] 
\ 3-P ^ 6-2P ^ ^ 


J 


niH = 


(i+iVe) J 



'92 



decreases 

i-p 



this range of q2, the value of {\6i 

as {q2e^°Y~^ in the massive case — or as {q2e"^}" " m 
the massless case — until the onset of large-(72 behavior 
at 52 = 1- The values of the exponents containing P 
can be seen from Eqs. ((50)) . ((56)) . and ([58]) . (Recall that 
P = 2 in the massless case, and the range of P for the 
massive case is shown in Fig. [T^l) The value q2 — ^ ef- 
fectively serves as a cut-off because of the rapid fall-off 
of ( I 6(j)^ I /iJ^ with increasing q2 in the large q2 re- 
gion. We find that the actual height of the major peak 
corresponding to the one normalized to unity in Fig. [TO] 
can be approximated as J e^'^"^)^'' in the massive case, 
and as J e*^'^^^-'^" in the massless case, with J having 
the same value, J ~ 0.01. Then we find that 



0.01 <mH< 1.49 
(01^0^10),, 



^infl 



dq2 J q2 



2-P 



•niH = 



(O|^0^|O)..c 



dq2 J ^2 



i-p 



(60) 



The resulting configuration space dispersion is given in 
TABLElH From this, one can deduce that the small mass 
limit, for which P ^ 3, gives the same result for \S(p\, as 
one has in the massless case, namely 



Hi, 



10 



(61) 



From the behavior discussed in the caption to Fig. 1131 one 
can show that this continuity also holds for each mode 
separately. 



X. SPECTRAL INDEX 

We say that a given mode k of the perturbation 
field 6(1) is crossing the Hubble radius at the time when 
k/{a{t)H{t)) = 1. For larger values of k/{a{t)H(t)) 
(shorter wavelengths) the mode is said to be inside the 
Hubble radius, and for smaller values it is said to be out- 
side the Hubble radius. Modes in the intermediate-(72 
range, as defined after Eq. p7p . exit during inflation to 
eventually re-enter the Hubble radius at some time after 
inflation has ended. Using our composite scale factor, we 
note that after a few e- folds of inflation, the quantum per- 
turbations that are exiting the Hubble radius are found 
numerically to satisfy — 



k^D{mi{) ' 



(62) 



where it will be recalled that 
given by Eqs. ^ and ^ 



as 



The variable Dirnn) is a 
constant of order 1 that we have evaluated numerically 
to be — 



D{mH 0) = 1.00, 
D{mH = Vol) = 1.04, 
D{mH = 1) = 1.45. 

One flnds approximately that Dimn) — (1 



(63) 
J)2. A 



horizon, the value of ji/'fe 
of 1 /2 that given in Eq 



few Hubble times after the mode has exited the Hubble 

approaches a constant value 
when TTiH — 0. In the case 
when ruH ^ 1, the value of {ipkl approaches about 1/2 
the value in Eq. after a few Hubble times, but then 
decreases very slowly over many Hubble times. 
The scalar spectral index rig is defined by 



n,, = 1 



d In Vs^ 
dink 



(64) 



At the time of exit from the horizon, or at several Hubble 
times after exit, we see that the spectrum defined by 
Eq. (|5T|) satisfies the proportionality 



Vs^ cx e 



(65) 



with 



IS 



oc k . Therefore the spectral index n 

1 for constant H as measured at any given number of 
Hubble times after horizon exit during inflation. This is 
true for the range of values oimn considered in Eq. ((63|) . 

The modes that exit the Hubble radius at early times 
before, or shortly after, the scale factor begins grow- 
ing exponentially are not described by Eq. (|62|) . These 
modes, which are primarily in the small-g2 region, are 
not scale- invariant. The gravitational perturbations in- 
duced by these modes could reenter the Hubble radius 
during the matter-dominated or dark-energy dominated 
stages of expansion, and would appear as long wavelength 
modes having a spectrum that is not scale invariant. 
If the total number of e-folds of inflation is sufficiently 
small, it would be possible to observe a breaking of the 
scale-invariance of the large-scale structure of the uni- 
verse at sufficiently large scales. See Figs. SI HI and [HI 
in which the breaking of scale-invariance is seen in the 
small-(72 region. The curves of \og{\(jq^ p) have a slope of 
—2 in the scale-invariant region of intermediate values of 
q2 and a slope of in the region of small 52 values. 

Because the svaaX\-q2 modes of large enough wavelength 
exit the Hubble radius before evolving away from the 
early-time conditions specified by Eq. (jTO]) , we would ex- 
pect a massless inflaton to generate a spectral index of 
rig = 3 in the small-(72 region. (This follows from the fact 
that the square of the amplitude in Eq. (fTO)) goes as k~^ 
for mass.) Similarly, we would expect a massive infla- 
ton to generate a spectral index of ris = 4 in the small-(72 
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region, if mjy satisfies q2 <C exp(— A^e) 'ttih in that re- 
gion. If scale-invariance continued indefinitely for large 
wavelength modes, the dispersion would be infrared di- 
vergent, so we expect an end to scale-invariance to occur 
at very large length scales for any reasonable spectrum 
of inflaton perturbations. Let us examine this breaking 
of scale invariance for the a{t) we have been considering. 
The modes responsible for galaxy-size structures today, 
left the Hubble radius approximately 45 e-folds before 
the end of inflation [H, pp. 285], so if were not too 
much larger than 45, we would expect it to be possible 
to observe the breaking of scale invariance at sufhciently 
large scales within our observable universe. 



It is of interest to understand the behavior of the 
coefficients E{k) and F{k) in the expression, given in 
Eq. ([SS]) . of the mode functions as a superposition of 
Hankel functions. These coefficients determine the quan- 
tum state of the inflaton perturbation field 54>. Because 
the modes evolve separately in this model, the state vec- 
tors can be considered for each mode separately. For the 
modes that exhibit a scale-invariant spectrum, we will 
see that the quantum state is the Bunch-Davies vacuum 
during inflation and that it evolves naturally from the 
initial Minkowski vacuum that we choose in the early 
time asymptotically flat segment of the scale factor. We 
find numerically that |-E(fc)| ~ 1 and F{k) ~ for modes 
of intermediate-q2, which are the modes that exit during 
the exponential expansion of our composite scale factor. 
Some sample values of E{k) and F{k) for our composite 
scale factor are given in TABLE Ull where the values of 
q2 listed assume a value of A^e = 60. The first row of 
data is in the small-g2 region. The second row of data is 
roughly at the interface between small- and intermediate- 
q2, and the fifth row of data is at the interface between 
intermediate- and large-q2- 



Thus, for intermediate-q2 modes (the ones that exit the 
horizon after inflation has started), and large-q2 modes 
(the ones that never exit the horizon) the value of |£'(fc)| 
is essentially 1 and the value of |-F(fc)| is essentially 0. 
This corresponds to the Bunch-Davies vacuum [23| in de 
Sitter spacetime. Since E{k) and F{k) are constant in 
time, this fact must be a consequence of our choice of 
state in the initial asymptotically flat spacetime at very 
early times. We took that state to be the Minkowski 
vacuum, having no particles present at early times. For 
sufficiently high frequency modes, there are no particles 
created by the initial expansion prior to the time ti when 
we join it to the inflationary segment (i.e., prior to the 
time when inflation begins in our model). The rate of 
expansion at this joining is equal to iJinfl, so one would 
expect few particles to be present at time ti in modes 
for which the momentum satisfies k/a(ti) > i/infi- This 
condition implies that the number of created particles 
present at the joining to inflation is negligible at values 
of 92 larger than those in the small-g2 region. Thus, it 



TABLE IL Sample Hankel Coefficients for Ne = 60 



52 




r (K) 




— 4H7S 26 — 5? 


— 4S7S 26 — .57 


^0-26 


-1.0915 - 0.017581i 


-0.18204 - 0.39819i 


^Q-23 


0.99924 - 0.038888i 


-1.68347 X lO""' + 4.37502i x 10"" 




-0.52204 + 0.85292j 


3.73431 X 10"" - 2.28562j x 10~" 


1 


-0.9082 - 0.41854i 


2.09333 X 10"^'^ - 2. 78876* x 10~^'^ 


10^ 


0.16978 - 0.98548i 


6.61571 X 10"*^ - 1.139902 x 10"*^ 



is reasonable^ that the initial Minkowski vacuum goes 
over into what is essentially the Bunch-Davies vacuum 
for values of k that satisfy this condition at the time 
inflation starts in our model. 

The large-(72 modes do not exit the Hubble horizon 
before inflation ends with the joining to the late-time 
asymptotically Minkowskian region at t2 in our model. 
Hence, for those modes one still has at late times that 
k/a(t2) > -ffinfi, so in those modes there are no created 
particles of the field 6(j), i.e., no perturbations of the in- 
flaton field are created. This behavior can be seen in 
Figs, m [51 andO where a rapid fall off in |/3fcp is seen 
for the large-92 region. By a similar argument to that 
in the previous paragraph, we would expect the number 
of particles of (50 to be negligible in the final Minkowski 
region for these modes k that did not exit the Hubble 
horizon during infiation. For these modes, the Bimch- 
Davies vacuum goes over into the Minkowski vacuum at 
late times. 

On the other hand, for intermediate-g2 modes 
a{t2)/a[ti) is sufficiently large that fc/(a(t2)^finfl) 1- 
Therefore, we would expect the intermediate-(72 modes 
to contain created particles of the 54) field in the final 
Minkowski space. This can be confirmed from Figs. [H [51 
andm 

For these intermediate values of 172, the value of the 
argument z of the Hankel functions is small in the late 
stage of inflation. Therefore, one can use the asymptotic 
form, 



r{v)Y fi 



(66) 



in Eq. ^ to find that IV^fc] ~ a"- 
^-3-,-2t)^ where z = k/{a{t)Hinfi) 



and I 



\J (9/4) — mj. From Eq. (f66l) . we see that |V'fc| is 
constant when mu — and decreases with time as 

a(i)(^ 



9_4m^2_3) > > 0. 

Using this small argument approximation with 



For a discussion of the Bunch-Davies vacuum in connection with 
the energy- momentum tensor see [25II . 
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Eq. (Uni) leads to 

din A: 

Then Eq. ([M)) gives a spectral index of 



9 - 4mJ, 



9 



(67) 



(68) 



When TJifj = 0, this gives the scale- invariant spectrum 
of inflaton fluctuations corresponding to = 1. The 
slow-roll approximation will reduce the value of Us be- 
cause i?infl decreases with time. This easily could be 
incorporated into our model, but will not be pursued 
further here. 

For 3/2 > > 0, Eq. ((68)l would seem to give > 1, 
but this is because in obtaining that equation we have 
used the same given t for all modes, which means that z 
is proportional to k for each mode. However, the pertur- 
bation spectrum of 54> is believed to induce the spectrum 
of scalar perturbations of the metric when the modes are 
at a given number of wavelengths outside the inflationary 
Hubble horizon [H, [2^ . This corresponds to the same 
value of z, not for each mode. When that is taken 
into account, the quantity dlnT^/lnfc = 0, thus giving 
the scale-invariant value Us = 1 for any value of mn- 

Because the gravitational perturbations are effectively 
massless, if one were to include them and propagate them 
to late-times using the scale-factor a{t) in our model, 
one would find the same scale-invariant spectrum as we 
obtained for an inflaton perturbation field with mu ^ 1. 

The graph of the inflaton spectra at late times for in- 
termediate q2 in Fig.[TT] has the spectral index of Eq. ([55]) 
because the joining to the asymptotically flat late time 
region occurs at a given time t2- For massless inflatons, 
as for massless metric perturbations, the fact that infla- 
tion ends approximately at a given time does not affect 
the scale-invariance of the spectrum.^ 



XI. DENSITY PERTURBATIONS 

In this section, we assume as usual that the inflaton 

perturbation fleld sets up scalar perturbations of 

the metric shortly after their wavelengths exit the Hubble 
horizon. We will consider two scenarios. 

In the first scenario, the modes of the quantized infla- 
ton field directly give rise to the scalar metric perturba- 
tions, without taking into account any regularization or 
renormalization of the ultraviolet divergences of the dis- 
persion (or variance) of the quantized inflaton perturba- 
tion fleld, beyond the Minkowsi vacuum energy subtrac- 
tion. If there are no further subtractions during the time 



when the scalar metric perturbations are induced, then 
the dispersion of the quantized inflaton field has an ultra- 
violet divergence. It is only in the late-time Minkowski 
space that this subtraction would give a finite dispersion 
for the quantized inflaton field. ^ Thus UV divergences 
are ignored in this first scenario, presumably because the 
modes relevant to observable scales today are assumed to 
be unaffected by renormalization of the ultraviolet diver- 
gences in the dispersion. However, it is not at all obvious 
that this is correct. This will become evident when we 
consider the second scenario. 

In the second scenario, the modes of the quantized in- 
fiaton field will again induce scalar metric perturbations, 
but we will regularize or renormalize the dispersion of the 
ultraviolet divergences of the quantized inflaton field, in- 
cluding the renormalization terms that are required for 
the dispersion to remain finite during the expansion of 
the universe, and not just the terms required in the late 
time Minkowski limit. We assume that at a time shortly 
after exit from the Hubble horizon, the scalar perturba- 
tions of the metric are induced and may be treated as 
essentially classical after that time [1, [l^]. It is most 
straightforward to use adiabatic regularization 0] to find 
the dispersion because it gives the result explicitly for 
each mode of the quantized inflaton perturbations. It 
should not make a significant difference to our result if 
one were to use, for example, Hadamard point-splitting 
regularization to obtain the dispersion ^j. One way to 
understand the rational for adiabatic regularization is 
analogous to that originally given in 0,0] in which a mea- 
suring instrument was considered. In the present case, 
one may suppose that the classical gravitational metric 
acts like a classical measuring instrument that cannot 
support the ultraviolet divergences that arise from the 
short-wavelength modes of the quantized inflaton fleld. 
The natural subtraction procedure that arises from this 
assumption is embodied by adiabatic regularization. 

We now evaluate the expected dispersion of the scalar 
metric perturbations by both methods for the simple in- 
flaton potential quadratic in the inflaton field. We work 
with the case when iJinfi is constant. As in 0, [H]) we 
find that when regularization is taken into account, the 
relevant metric perturbations induced by a given value of 
-ffinfi are considerably smaller than they would be if no 
regularization of the quantized inflaton fleld were taken 
into account. (By "relevant modes," we mean those re- 
sponsible for temperature variations in the CMB and 
large scale structure that we observe today.) Put an- 
other way, when regularization is taken into account, the 
value of i?infl responsible for the these observed variations 
must be larger than it would be without regularization. 

The basic reason for the smaller dispersion found in 



The scale-invariance of the initial spectrum of scalar perturba- 
tions of the metric when inflation ends, thus places an upper 
limit on the graviton mass, at least in principle. 



This can be seen, for example, by looking at the terms in the 
adiabatic subtraction that vanish when a and a are — those 
terms are divergent when integrated over all modes. 
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the second scenario for the scalar metric perturbations is 
that the metric perturbations are induced, and acquire 
classical properties, within about a Hubble time after the 
inflaton perturbations have left the horizon. The disper- 
sion of the classical scalar metric perturbations depends 
on the value of the dispersion of the inflaton field at the 
time when the perturbations are induced. Because that 
time is soon after the inflaton perturbations have exited 
the Hubble horizon, the adiabatic subtractions are quite 
significant. After that time the induced scalar metric 
perturbations are treated as classical, requiring no fur- 
ther subtractions. 

If the scalar metric perturbations were induced only at 
the end of the inflationary era, then the relevant modes of 
the inflaton perturbation field would have left the Hori- 
zon a large number of Hubble times earlier and the adi- 
abatic subtractions would be small. In that case, the 
adiabatic subtractions would give only a small difference 
from the first scenario. However, in the second scenario 
we are assuming the metric perturbations are induced 
soon after horizon exit. 

In the case when mu^ <C 1, we find numerically that 



the duration and manner in which inflation ends does not 



have much affect on the value of 



after a given 

mode in the relevant range of wavelengths has exited the 

2 



at late times 



Hubble radius. Thus, the value of 
is approximately equal to the value it has soon after leav- 
ing the Hubble horizon. At late times, Eqs. ([T7]), 
and (|22p show that the expectation value is 



^(un) 



1 



1 



'2/'^2/ 



(|afe|' + |/?fe 
(l + 2 



(69) 



The value of 



obtained from Eqs. ((T6|), (fTT]) . 

and (|62p. however, is un-renormalized. In the late time 
Minkowski spacetime, the renormalization consists of 
subtracting the vacuum contribution corresponding to 
the late time Minkowski spacetime for which the scale 
factor has the value 02/. This is equivalent to replacing 

(1 + 2 l/Jgal^) by 2 l/Jgal^ in Eq. ([69]) . as was done to obtain 
Eq. (gll). Thus, at late times 



(70) 



For intermediate values of qi in the massless case, it 
follows from Eq. PO)) that at late times 



1 



792 



(71) 



We write the renormalized spectrum of S(j){x) in terms 
of the modes of the renormalized expectation value of 
the dispersion (or variance), (0 |((50(x))2| 0)ro, which are 



defined in Eq. (p4)l . This gives the spectrum as 

3 



27r J 



(72) 



From Eqs. ^ and we obtain, in the continuum 

limit, the late-time renormalized inflaton perturbation 
spectrum 



87r2 



(73) 



From our discussion of Fig. [131 this result also holds in 
the massive case when m//^ <^ 1. Here, in the range 
where m? ^ fc^/a2/^, we have l/^gal^ °<: (72 As can 
be seen in Fig. 1131 this region is continuous with the 
massless case given by Eq. (|7ip for fc^/a2/^ ^ rri^. This 



continuity tells us that the region of \(3q^\ 
to q2~^ must follow 



and thus 



ruH _3 
—92 , 



proportional 



(74) 



87r2 



(75) 



for both the massless case and the case where niH^ <C 1. 
This is the result for the spectrum in the first scenario, 
in which only the late-time Minkowski vacuum subtrac- 
tions are considered and the renormalization of the in- 
flaton dispersion during the expansion of the universe is 
ignored. 

However, by late-times it is usually assumed that the 
inflaton perturbations have already induced the scalar 
metric perturbations that will give rise to the initial scale- 
invariant plasma oscillations in the reheated universe. 
Because the universe is still expanding rapidly at the time 
that the relevant scalar metric perturbations are induced, 
one should consider the effect of renormalization of the 
inflaton perturbations in the rapidly expanding universe 
on the amplitude of the metric perturbations that they 
induce. Following 0, [3j we take the scalar metric 
perturbations to be classical in nature shortly after they 
are induced. Therefore, in our second scenario, it is the 
renormalized spectrum of the quantized inflaton pertur- 
bation field shortly after their wavelengths exit the Hub- 
ble horizon during inflation that we use to determine the 
magnitude of the spectrum of the induced scalar metric 
perturbations. The latter propagate as classical quanti- 
ties during the long remaining duration of inflation, and 
therefore their amplitudes remain almost constant up to 
the time that they induce acoustic plasma oscillations 
after reheating. The adiabatic subtractions coming from 
renormalization are significant at the time that the scalar 
perturbations of the metric are induced, and hence have 
an important influence on the magnitude of the metric 
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perturbations and consequently on the CMB anisotropies 
and the large scale structure observed today. Slow-roll 
inflation will alter the scale-invariant spectrum, as con- 
sidered with the effect of renormalization in Q. The 
effect of renormalization is found there to be important 
for interpreting the existing measurements of the CMB 
anisotropies. 

Next, we consider the effect of renormalization on the 
amplitude of the spectrum of scalar metric perturbations 
induced by the quantized inflaton field during inflation. 
We give the results for both scenarios; the first, in which 
only the Minkowski subtractions are made and the late- 
time dispersion in the final Minkowski space is used; and 
the second, in which the renormalization is carried out 
near the time when the wavelength of the mode first ex- 
ceeds the Hubble radius during inflation. 

The curvature perturbations (or scalar perturbations) 
of the metric [1^ are defined by 



(76) 



where (j) represents t/)*^"^ in the present section. In [2^ . 
the possible effects of renormalization by adiabatic regu- 
larization or point-splitting were not considered, and the 



unrenormalized value 



was used. However, here 



we will use the renormalized value, as it is clear that at 
late times one must subtract the Minkowski vacuum con- 
tribution. This subtraction is already included in adia- 
batic regularization (and in point-splitting) when applied 
to Minkowski space. However, when d or d are nonzero, 
the simple Minkowski space subtraction is not enough 
to yield a finite value for the dispersion of the inflaton 
fluctuation field in the expanding universe. As noted in 
@i 0, Q I the effect of renormalization in the expanding 
universe may have a significant effect on the magnitude 
of the inflaton dispersion spectrum. As discussed in the 
previous paragraph, we take the relevant time at which 
the magnitude of the inflaton perturbation fleld induces 
scalar metric perturbations to be shortly after horizon 
crossing. Therefore, in the second scenario we use the 
horizon crossing time as a characteristic time at which to 
evaluate the renormalized spectrum in Eq. (|76p . One will 
obtain a similar result as long as the characteristic time 
is taken within a few Hubble times after horizon crossing. 

The spectrum, Vn = {L/2Tr)^4iTk^7ik^ , of scalar met- 
ric perturbations in the continuum limit now takes the 
form 



(77) 



The difference between the two scenarios boils down to 
the time at which the renormalized dispersion is calcu- 
lated. 

So far, our parameterized scale factor has not been 
linked to any particular potential or model of inflation. In 
what comes next, we choose a simple potential that was 



found to be in good agreement with the 3-year WMAP 
data llj. For our example, we use the quadratic chaotic- 
inflation potential [2^ 



2 ^ 

The two slow roll conditions are [l| 

SttG, 



-V, 



and 



dV/d(f> 
3H ■ 



(78) 



(79) 



(80) 



We combine these two slow roll equations with the 
present quadratic potential to find 



1 

3 ^/8^■ 



(81) 



Then with the quadratic potential, in the first scenario, 
we find from Eqs. ([75]) and ([77|l with only the late-time 
Minkowski subtraction, that 



ifi 



2m2 [l/Vs^G 
which can be written as 

3 / -ffinfl 



87r2 ' 



(82) 



Vn = 



167r2m//2 \ 2.436 X lO^^GeV 



(83) 



In the second scenario, we take account of the further 
subtractions necessary to renormalize the inflaton dis- 
persion when the time-derivatives of a(t) cannot be ne- 
glected and obtain for the inflaton dispersion spectrum. 
We use the result given in Q for the renormalized spec- 
trum, Vg^iTe) . As noted before, this should be evaluated 
shortly after the mode has exited the Hubble radius. For 
simplicity, we use the time of exit to characterize the 
order of magnitude one will obtain. This gives 



32 7r2 



47r 



H 



(1) 



■(1) 



iff 



H 



16 



from which one obtains, 
3 



Vn 



^7/2 



& 



(84) 



64 7r2m/f2 \^ 2.436 x W^GeV 

2 



X 47r 



H 



(1) 



■(1) 



8mH^' + iiniH^ + AGmH^ + 16 



(rn/f2 



1 



,7/2 



(85) 



18 



After inflation has ended, at the time when once again 
k/{a{t)H{t)) — 1, the curvature perturbations reenter 
the Hubble sphere. At that reentry time, their amph- 
tudes are related to the amplitudes, Spk, of the density 
perturbations by relations of the form [23|, 



Spk 
P 



(86) 



where x is 2/5 if reentry occurs during the matter- 
dominated stage, and is 4/9 if reentry occurs during the 
radiation-dominated stage. We use the value 2/5 in the 
rest of this section. Then, in the first scenario, the spec- 
turm of density perturbations created by the scalar met- 
ric perturbations at the time of reentry is 



5 



100 TT^ m/f" 



2.436 X 10^»GeV 



(87) 



In the second scenario, assuming that the magnitude of 
the scalar metric perturbations does not change much 
from the value it has at the time that the metric pertur- 
bations are induced to the time of reentry, we have 



Vs 



400 7r2TO^2 



X i-K 



H 



(1) 



2.436 X lO^^GeV 

2 

(1)' 
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See TABLE IIIII for sample values of the density con- 
trast, Vs, given by Eq. ([57]) in the first scenario, and see 
TABLE lTVl for sample values of the density contrast given 
by Eq. ([88]) in the second scenario. In both cases, the po- 
tential is the quadratic one of Eq. ([75]) . The quantity Vs 
in our table is related to the measured quantity denoted 
by in the WMAP five-year results 30,, see their Ta- 
ble 2]. Their measurements give for the most likely value, 
= 2.41 X 10~^. In our present paper, A^ is denoted 
by V-R,, and from the first line of Eq. (87), we obtain its 
observed value: Vs = {2/bfVn ~ 3.86 x 10-^°. We have 
chosen H to give values of Vs in the first row of each 
table that are near the observed value of Vs- 

A comparison of TABLE |TV] with TABLE Hn] reveals 
that adiabatic regularization leads to a constant non- 
zero value of Ps for small values of mn- Thus, for the 
quadratic potential, all the values of mn shown for adi- 
abatic regularization in TABLE IIVI give a value of Ps 
that agrees with the observed value when H w 10^^ 
GeV. In TABLE IIIII we see that Ps grows to as niH 
decreases. For the quadratic potential it approaches in- 
finity in the limit as mn — >■ without adiabatic regular- 
ization, but approaches a constant non-zero value when 
the full adiabatic regularization is used during the ex- 
pansion. This can be confirmed by expanding the right- 
hand-side of Eq. in powers of tuh and taking the 



TABLE III: Late- Time Method: Vs for V = im^(^^ 





H = lO^GeV 


H = 10" GeV 


H = lOi^GeV 


niH = 0.0001 


5.12 x 10"^^ 


5.12 X 10-1° 


5.12 X 10-* 


niH = 0.01 


5.12 X 10^^^ 


5.12 X 10'" 


5.12 X 10-1^ 


rriH = 0.1 


5.12 X 10"^** 


5.12 X 10-1® 


5.12 X lO-i'' 


rriH = 0.25 


8.20 X lO"^'' 


8.20 X 10-" 


8.20 X lO-i'"' 


TABLE IV: Adiabatic Regularization: Vs for V = ^m'^cj)^ 




H = lO^-'GeV 


H = lOi'^GeV 


H = lOi^GeV 


niH = 0.0001 


4.60 X 10"^^ 


4.60 X 10-1° 


4.60 X 10-** 


niH = 0.01 


4.60 X 10"^^ 


4.60 X 10-1° 


4.60 X 10-* 


niH = 0.1 


4.64 X 10"^^ 


4.64 X 10-1° 


4.64 X 10-* 


niH = 0.25 


4.79 X 10"^^ 


4.79 X 10-1° 


4.79 X 10-* 



limit as mH 0. The fact that the value of Ps is almost 
constant in TABLE HVl for all the values of mn < 0.25 in 
each column, implies that one could predict the value of 
i?infl for this, and perhaps also for other inflaton poten- 
tials, when adiabatic regularization is taken into account. 



XII. REHEATING 

When we maintain continuity of a(t), d{t), and d{t), 
the particle number in a given mode is proportional to 
q2~^ for large enough values of (72 • [See Eq. (|42l) .] 

When the end of inflation is gradual, we found (as seen 
in Fig. that beyond 52 — 1, the quantity |/3g2p f^l^^ 
off as q2~^- For small mn, using the behavior of l/Sgal^ 
summarized in Eqs. (HOI), and for the case of 

a gradual end to inflation with a2//ai/ > 2, we find (by 
blue shifting the temperature at late times back in time) 
that at the time when the inflationary expansion is joined 
to the final asymptotically flat segment of the expansion, 
the energy-density has an effective temperature that is 
slightly larger than the Gibbons-Hawking temperature 
|3ll | of i?infl/27r. Particle creation by the expanding uni- 
verse was also investigated in [l^ as a possible cause of 
reheating, with a similar result to the present one we find 
for a gradual end to inflation. 

We next use our results to investigate the effective 
temperature if there is a fairly abrupt end to infla- 
tion. With a fairly abrupt end to inflation, for which 
0-2/ /aif ~ 1, the range in (72 where j/Sg^j oc q2~^ 
only begins until a large value of g2, which we denote 
by 92cut-off- From our numerical results, we find that 
<Z2cut-off' — 0-2/ / (0-2/ — 0,1 f)- We define the region between 



1 < 



92 



< 



92cut- 



"stretched" region. In this 

-2 



(72 As an example, see 



-off as the 
"stretched" region, IPqj'^ cx 

Fig. [51 where Q2cut-off — 10** and a2//(a2/ - ai/) ~ 10''. 
The "stretched" region is a result of particle production 
caused by the rapid fall-off in H{t) as inflation ends. 
With sufficient stretching, the particle number per 
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mode, |/3g2 1 , in the "stretched" region is proportional 
to 92"^ f'-'^ value of rriH < 3/2. In the intermediate- 
q2 region, the behavior of j/Sg^ |^ is governed by the power 
P defined in Eq. ([55]) and shown in Fig. [T^l As is evident 
from Fig. [HI the stretched region extends over a much 
larger range of values of (72 than do the intermediate- 
and small-q2 regions. Because the a(t) we defined is only 
across the joining points, 1/3,2 1^ off as 92"^ for 
<Z2 > 'Z2cut-off- For an a{t) that is C°° for all t, the rate 
of fall off for q2 > (72cut-off would be much faster. There- 
fore, in our calculation below of the total energy density, 
we neglect the contribution to the energy density from 
this UV-range of (72- (However, in the case, one can 
show that the contribution of this UV-range to the total 
energy-density is of the same magnitude as the contribu- 
tion of the stretched region; so in that case the energy 
density would be twice the value we obtain below.) 

Therefore, for a fairly abrupt transition with 
o-2f/o-if — Ij the total energy density is 



1 /■a2/-Wi„fl g2cut-off I, 



E 

V / ~ (2TXaisY Ja2fH,^n ^2/ 

■^—521^^4(1 
1 27r2 

<j2cut-off „ 717-4 



'-dq2 



87r2 



dq2 



(89) 



In the present case, q'2cut-off — (02// [02/ - ai/]) > 1, so 
we find 



E \ ^in 



V 



( V 

fi ya2f-aif J 



167r2 



(90) 



It follows that a fairly abrupt end to inflation can lead to 
energy densities that correspond to an effective temper- 
ature large with respect Gibbons-Hawking temperature. 
This would certainly be large enough to reheat the uni- 
verse if some process caused a fairly rapid end to infla- 
tion. A more gradual end to inflation could also produce 
a temperature sufficient for reheating if the value of i/infl 
were large enough, as pointed out in [l^ . 



XIII. CONCLUSIONS 

We have set up an exactly solvable model of exponen- 
tial inflation joined to initial and final asymptotically flat 
expansions of the universe. There are eight adjustable 
parameters, three of which pertain to the initial asymp- 
totically flat segment, three to the final asymptotically 
flat segment, and two to the inflationary stage of the 
expansion. In the case when the scale factor a{t) is re- 
quired to be continuous with continuous first and second 
derivatives at the joining times, there are three indepen- 
dent conditions on the parameters, leaving five adjustable 
parameters. Among these are the rate of inflation, -ffinfl, 
the number of e- foldings of inflation, N^, as well as three 



parameters characterizing the asymptotic regions. The 
parameter iJinfl is constant in this model, and N^, can 
be arbitrarily large. Our principal focus in this work has 
been to use this exactly solvable model to look for generic 
behaviors that have not been previously studied. In this 
way, we have found a number of new effects that may 
be relevant to the inflationary universe. We summarize 
these as follows. 

We find that if the quantized inflaton fluctuation field 
is in its vacuum state in the early time asymptotic 
Minkowski space, then the inflaton fluctations (which 
correspond to particles of a quantized scalar field) are 
created by the expanding universe in the manner thor- 
oughly analyzed in [6, 7]. We show that the creation of 
these particles, or infiaton fluctuations, leads to a scale- 
invariant spectrum^ of created particles, or fluctuations, 
for the range of momenta or wavelengths that are relevant 
to the present universe when Ne is sufficiently large, for 
example, larger than about 60. For an initial Minkowski 
vacuum, we find that for modes having wavelengths that 
are not so large that they leave the Hubble sphere within 
the first few e-foldings of inflation, the Minkowski vac- 
uum state evolves into the Bunch-Davies de Sitter vac- 
uum after a few e-folds of infiation. So the results for 
such modes of the inflaton fluctuation field should agree 
with those obtained by assuming the de Sitter vacuum, 
as is usually done in treatments of inflation. 

If the total number of e-foldings of inflation prior to 
the exit of all the observationally relevant modes is not 
sufficiently long, then our results would imply that the 
longest wavelength perturbations that may come within 
the range of observation should have a spectrum that 
deviates in a certain way from scale invariance. This de- 
viation is a result of the fact that the Minkowski vacuum 
has not had time to evolve to the de Sitter vacuum before 
the modes exit from the Hubble sphere. 

We find that if there were perturbations (i.e., particles) 
present at early times, then their effects can be propa- 
gated all the way to the end of inflation, by means of 
"stimulated emission" from the vacuum. These pertur- 
bations present prior to the start of inflation would then 
have observable effects. This is an important new effect 
that we are considering in a paper{32] in preparation. 

We investigated the spectrum of scalar perturbations 
of the metric that would be induced by the inflaton fluc- 
tuations, and we found the spectrum of initial density 
perturbation that these would produce after the end of 
inflation. We found a signiflcant difference between the 
properties of the density perturbations that one would 
obtaiin by taking account (a) only of the vacuum sub- 
traction in the flnal Minkowski space, as compared with 
(b) subtracting the full set of adiabatic regularization 
terms at a time within a few e-foldings after a mode cx- 



* We are not incorporating the slow-roll parameters, which of 
course would modify the scale-invariance. 
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its the Hubble sphere during inflation, as in Q. A new 
result wc find is that in case (b), the density perturba- 
tions have a non-zero value that is nearly independent 
of the inflaton mass for masses < 0.25.^. This type 
of behavior may permit one to use our knowledge of the 
initial density perturbations to make a prediction about 
the value of H during inflation. This mass-independence 
does not appear if one only uses the Minkowski space 
vacuum subtraction. 

We also considered the energy density that would be 
created by the change of the scale factor a{t) at the end of 
inflation. We found agreement with earlier treatments for 
a gradual end to inflation. We also were able to estimate 
the energy density created if there were a fairly abrupt, 
but smooth, end to inflation. The effective reheating 
temperature produced could be very high. We give the 
value of the created energy density as a function of a 
parameter characterizing the abruptness of the end to 
inflation. 
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APPENDIX A: BOUNDARY MATCHING 
CALCULATION 

Given the values of ipki and -0^]^, and the matching 
conditions 



Ahiih) + Bh2{ti) 
Ah[{ti) + Bh'^iti) 
Cgi{t2) + Dg2{t2) 
Cg[{h) + Dg'^ih) 



V'fe(ii) = 0fei, 

Ahi{t2)+Bh2{h), 
Ah[it2)+Bh'2it2), 



(Al) 



we wish to calculate the constant coefficients C and D in 
terms of the functions hi{t), h2{t), gi(t), and g2{t)', and 
the values of i/'fei, V'feii ^i; and t2. (Here prime denotes 
derivative with respect to t.) Rearranging the first two 
matching conditions leads to 



B -- 

A- 



V'fei - Ahi 



(A2) 



kl 



BK 



t=tl 



^ At least for the quadratic inflaton potential that we considered, 
but probably also more generally 



Combining these two equations leads to 



A 



B = 



h[h2 — hih'2 

i^'kihi - ipkiK 

h'2hi — h2h[ 
At the time, ^2, we have: 

Mt2) = Ahi{t2) + Bh2{t2) 
i^'kih2 - TpklK 



t=ti 



(A3) 



and 



h'ih2 — hih'2 

/l2^1 — /l2^'i 



ijiit2) = Ah[{t2) + Bh'2it2) 
V'fel^2 ~ ^Pklh'2 



hl{t2) 



h2{t2) 



(A4) 



+ 



h'^h2 — hih'2 
h'^hi — h2h[ 



t=tl 



h'2{t2) 



t=ti 



(A5) 



Let us also define ipk2 = '4'k{t2) and = V'fcl^z)- In 
terms of ipk2 and the last two boundary conditions 
in Eq. (|A1[) become 



C 



D = 



i^'k292 - 1pk2g'2 

g'ig2 - gig'2 
^'k29i - i^k2g'i 



5251 



525i 



(A6) 



Substituting for ^pk2 and -0^2 yields 

[Ah[ + Bh'2]g2 - [Ahi + Bh2]g'2 



D = 



5i52 - 5152 
[Ah[ + Bh'2]gi - [Ahi + Bh2]g[ 
.9251 - .9251 



(A7) 



Finally, expressing A and B in terms of the given values 
of 0/ci and ■0^j^ specified at ti leads to 



C = 



1 



(5l52 - 5i52)t=t2 

■0fei^2 - i'kiK 



h[h2 



^'kihi 



— hih'2 

i>kiK' 



{h[g2 - /ii52)t 



t=ti 



/ij/ii — h2h[ 



(^252 - /i252)t=t2 



t=tl 



(A8) 



and 
D 



(5251 - 929'i)t^t2 
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h[h2 — hih'2 

-0^1 /ii - i^kiK 



■ti 



(A9) 



which are the combined joining conditions for and t/i^. 
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